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ABSTRACT 

Following the discovery of a new radio component right before the GeV 7-ray detection since 
2008 August by Fermi Gamma-ray Space Telescope, we present a detailed study of the kinematics 
and lightcurve on the central sub-pc scale of 3C 84 using the archival VLBA 43-GHz data covering 
the period between 2002 January to 2008 November. We find that the new component "C3", 
previously reported by the observations with the VLBI Exploration of Radio Astrometry (VERA) , 
was already formed in 2003. The flux density of C3 increases moderately until 2008, and then 
it becomes brighter rapidly after 2008. The radio core, CI, also shows a similar trend. The 
apparent speed of C3 with reference to the core CI shows moderate acceleration from 0.10c to 
0.47c between 2003 November to 2008 November, but is still sub-relativistic. We further try to 
fit the observed broadband spectrum by the one-zone synchrotron self-Compton (SSC) model 
using the measured apparent speed of C3. The fit can reproduce the observed 7-ray emission, 
but does not agree with the observed radio spectral index between 22 and 43 GHz. 
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The radio source 3C 84 is associated with the 
giant elliptical galaxy NGC 1275 (z = 0.0176). It 
is well known that 3C 84 has a pair of compact 
radio lobe structures on the central 10-pc scale 
(IVermeulen et al.[ ll994l IWalker et aDll994l 120001 
Asada et al.l 2006 ) . Because of its proximity, it al- 
lows us to study the region within the central sub- 
pc scale using a Very Long Baseline Interferometer 
(VLBI). It is an ideal laboratory to investigate the 
formation mechanism of relativistic jet ultimately 
powered by super-massive black holes and the in- 
teraction between the jets and ambient matter in 
the galactic central regions. 

Thanks to the recent observations with Fermi/h AT, 
we have a new opportunity to explore GeV 7-ray 
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produc tion mechanism in misaligned radio-loud 
AGNs (lAbdo et al.ll20inh . With CGRO/EGRET, 
a few extragalactic radio gala xies, suc h as Centau- 
rusA ( Sreekumar et aljl999h. 3C 111 (lHartman et~al 



sponds to 0.24c - 0.55c, i.e. slower th an the speed 



2008h . and NGC 6251 dMukheriee et al.ll2002h had 



been alread y detected. Howev er, the detection of 
NGC 1275 (|Abdo et all l2009h by Fermi/LAT in 
2008 August is particularly noteworthy because 
NGC 1275 was not detected by CGRO/EGRET 
( Reimer et al.l l2003h . The flux density detected 
by Fermi/LAT is about 7 times higher than the 
upper limit of EGRET sensitivity. Intriguingly, 
the radio monitoring also shows the flux increase 
starting in 2005. The time variation of 7-ray flux 
density shows a similar trend with the radio flux 
density on the timescale of decades, implying the 
possible connection between the 7-ray emission 
and variable radio component. 

Overall spectral energy distribution (SED) of 
NGC 1275 from radio to 7-ray can be explained 
by the synchrotron-self Comp ton (SSC) mode l 
and the deceleration jet model ( Abdo et al. 20091 ) 
adopting sub-pc for the size of emitting region. 
The authors derived Lorentz factors of the emit- 
ting region F = 1.8 (0 = 25 deg, 5 = 2.3) for the 
case of SSC model and T varying from 10 to 2 
(0 = 20 deg, S varying 1.6 to 2.7) for the deceler- 
ating model, where and S are the jet angle to the 
line of sight and the beaming factor of the emit- 
ting region, respectively. Here we should stress 
that the size of the 7-ray emitting region adopted 
in the SED models is comparable to a VLBI com- 
ponent size. Therefore, it is essential to test the 
scenario of co-spatiality of GeV 7-ray and radio 
emitting regions by VLBI observations. 

In order to find the radio counterpart of the 
GeV 7-ray emitting region, we have conducted 
VLBI observations using the VLBI Exploration of 
Radio Astrometry (VERA) at 22 G Hz between 
2006 June 14 and 2009 April 24 (jNagai et al. 
(<2010h . lereafter Paper I). Surprisingly, we found 
that the monotonic increase of radio flux den- 
sity mainly originated in the newly born bright 
component C3 (Paper I). The measured projected 
speed of C3, however, was (0.23 ± 0.01)c on aver- 
age between 2007 October 24 and 2009 April 24 
in the central sub-pc to pc scale. Given previ- 
ously e stimated jet angle to the line of sight (11 ° 



55°: lAsadaet al.1 (j2006l ); iLister et alJ (|2009f ) 



I Abdo et al.l ( 2009 )). the de-projected speed corre 



of jet derived from the SED modeling (jAbdo et al 
20091 ). The result of VERA observation implies 
S < 2, and there seems to be a discrepancy be- 
tween the SED models and VLBI observations. 

A possible idea to explain this discrepancy is 
that a jet component with relativistic speed other 
than C3, located much closer to the central black 
hole (within the central sub-pc) , is responsible for 
the GeV 7-ray emission. If such a relativistic mo- 
tion in the jet of 3C 84 occurs only in vicinity of 
the black hole, the speed of C3 could have been 
relativistic soon after the ejection from the core, 
and then it underw ent to a rapid deceleration. 
Liuzzo et al. ( 2010l ) also suggested that jets arc 
relativistic in bright cluster galaxies but they de- 
celerate very soon (sub-parsec scale) because of a 
strong interaction with the ISM. In Paper I, we 
argued possible detection of the relativistic speed 
at the early stage of emergence of C3. However, 
this was not conclusive because of (1) the lack of 
spatial resolution of VERA at 22 GHz and (2) pos- 
sible screening by free-free absorption. 

In order to overcome the above potential prob- 
lems, higher resolution observations at higher fre- 
quency are required. We explore the possible rel- 
ativistic flow of C3 at the early stage of its emer- 
gence from the core with the archival VLBA data 
at 43 GHz. We selected only data obtained from 
2002 to 2008, i.e. considering earlier epochs than 
those discussed in Paper I. We also discuss the 
SED fit using one-zone SSC model involving the 
actual measured apparent motion of C3 and pos- 
sible ideas for incorporating the VLBI observation 
with GeV 7-ray emission. 

Throughout this paper, we adopt the follow- 
ing cosmological parameters; Hq = 71 km sec -1 
Mpc -1 , fl M = 0.27, and tt A = 0.73 (1 mas = 
0.353 pc, and 0.1 mas yr" 1 = 0.113c). 

2. OBSERVATION AND DATA ANALY- 
SIS 

2.1. The VLBA at 43 GHz 

First, we examined all of the (available) archival 
VLBA data at 43 GHz obtained between 2002 to 
2008 from the data archive system of the National 
Radio Astronomy Observatory (NRAO). Data re- 
duction was performed using Astronomical Im- 
age Processing System (AIPS). A priori ampli- 
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tudc calibration was done based on the measure- 
ments of the system noise temperature (T sys ) dur- 
ing the observations and the elevation dependent 
antenna gain provided by each station. In this 
process, we also applied opacity correction due to 
the atmospheric attenuation, assuming that the 
time variation of the opacity is not significant 
during each observation. Fringe fitting was per- 
formed with the AIPS task FRING. After apply- 
ing the delay and rate solutions, we combined each 
set of data into one channel across the 8-MHz 
bandwidth. Imaging and self-calibration process 
was performed using Difmap software package 
( Shepherd et al.lfl994j ) . The final images were pro- 



duced after iterations of CLEAN, phase and am- 
plitude self-calibration processes on successively 
shortened calibration time. Second, we excluded 
the entire datasets with relatively poor dynamic 
range (because of bad weather condition or lack 
of baselines) . Finally, we ended up with 28-cpoch 
datasets (shown in Tabic [T]). In almost all the 
archival data we analyzed, 3C 84 was observed as 
a calibrator and the typical observing time is less 
than 10 min. 

2.2. The Metsahovi at 37 GHz 

The single dish monitoring data at 37 GHz is 
adopt ed from Metsahovi quasar monitoring pro 



gram (jTerasranta et alj|1998t iLeon-Tavares et al 
201 if ). The observations were carried out with 
Metsahovi 14-m Radio Telescope. In order to com- 
pare the VLBI and single dish lightcurves, we se- 
lected the data from 2001 December to 2008 De- 
cember. Details of these observations and calibra- 



tion are described in iTerasranta et al.l (|1998l ). 



3. MODELING THE VLBA DATA 
3.1. Gaussian Model Fittings 

To obtain the flux density of each component of 
3C 84 jet, we modeled the source structure using 
the Difmap task modclfit. In the modelfit pro- 
cess, we fitted multiple elliptical Gaussian model 
components to the visibility data of each epoch. 
We adopted the fitting model components that 
gave a good fit to the data as judged by relative 
X 2 -squared statistics. In some cases, we adopted 
circular Gaussian models or point sources instead 
of the elliptical Gaussian if either of them were 
capable of reproducing a better fit. We consider 



the fit acceptable if the model components con- 
tain more than 95% of the total CLEANed flux. 
Resultant images mainly consist of three major 
components (CI, C2, and C3; shown in §4), al- 
though C2 is not identified in some epochs because 
of low signal-to-noise ratio. In addition to the ma- 
jor components, some faint minor components are 
fitted. The choice of adding a few minor compo- 
nents to the model is somewhat arbitrary, and it 
is difficult to maintain the consistency across the 
epochs. We note, however, that the choice of mi- 
nor components does not appreciably affect the 
fitting results of the major components. 

In order to check that the fit is not local min- 
imum, we confirmed that the result is unchanged 
significantly by choosing different starting Gaus- 
sian models. 

3.2. Positional Error of VLBA Images 

In general, it is difficult to measure a positional 
error of each component quantitatively from sin- 
gle epoch interferometric data. In particular, a 
systematic error, which arises from the deconvo- 
lution error, the presence of other nearby jet fea- 
tures and so on, is difficult to estimate separately. 
Hence the sampled interferometric-visibility on 
the (u, w)-planc is poor, and therefore systematic 
error of component position may be large com- 
pared to long-track observations. For these rea- 
sons, we estimated the component positional er- 
rors by examining the scatters in the positions of 
components (defined by the AIPS task MAXFIT) 
in images between two close epochs (within 10-day 
separation), such that the jet has approximately 
the same structure in both epochs. The appar- 
ent motion of lc will show 2.4 /zas/day positional 
change (< 0.01 of typical VLBA 43-GHz beam), so 
the motion of the component is negligible within 
10 days in this study. We have done this analysis 
on 11 pairs of images (Tabled]). 

Figure [T] shows the differences of position of C3 
for these 11 pairs. Each data point is normalized 
with the beam size averaged between the pair. The 
standard deviations along right ascension (af A ) 
and declination (c° EC ) normalized by beamsize 
(dheam) are .150 and 0.098, respectively. As in 



Barnes! (|l994h . the 100a% confidence interval of 
the standard deviation of population for statisti- 
cal ensemble i, at, is estimated from N sample of 
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statistical ensemble i as 



5&-i(«/2) p X 2 N -i(l-a/2)' v ' 

where X%( a ) 1S X 2 statistics for d.o.f. = N on 
which the event % 2 > x^- occurs with probability 
a. Wc apply this estimator for right ascension 
ensemble (i = RA) and declination ensemble (i = 
DEC). Given N = 21 and a = 0.95, the standard 
deviation of population for both right ascension 
(cr^ A ) and declination (o-p EC ) is 

0.106 < cr« A < 0.255, (2) 
0.070 < cr° EC < 0.166. (3) 

Hereafter, the positional error of each component 
is set as 0.255 #beam for right ascension and 0.166 
#beam for declination, where 0^ cam is the beam size 
on each epoch. 

4. RESULTS 

4.1. Overall Structures 

Figure [5] shows the total intensity images at 
43 GHz. Information on all images is also sum- 
marized in Table [3l All images are convolved with 
the circular restoring beam with a diameter of 
0.3 mas. After 2003 November 20, a new com- 
ponent C3 seems to emerge south of CI. This 
new component gradually becomes brighter with 
increasing separation from component CI. Emer- 
gence of C3 is also seen in 22-GHz observations 
(Paper I) around 2007, but the VLBA 43-GHz 
observations find that the component has already 
emerged in earlier epoch, because of thin opacity 
and/or higher resolution than the VERA 22 GHz 
observations. After the emergence, C3 is advanc- 
ing toward position angle ~ 170° measured from 
CI. Its direction differs from that of C2 with re- 
spect to CI by 40-50 degrees. No counter jet com- 
ponent is detected with a level of 2a throughout 
all epochs. Physical parameters of all fitted com- 
ponents are listed in Tables 21 \E\ and |6] 

4.2. Kinematics 

We use the AIPS task MAXFIT to define the 
(relative) positions of components at each epoch. 
Since absolute position in each image is lost at 
the self calibration and fringe fitting processes, we 



need to define the reference position for the ar- 
gument of kinematics. In this work, we will eval- 
uate kinematics with reference to the optically- 
thick radio core. Thus we attempt to measure 
the spectral index of each component to define 
the reference position. We measured the spec- 
tral index a for each component at 2008 Novem- 
ber 27, when 22 GHz and 43 GHz observations 
were carried out only one day apart. We define 
the spectral index between 43 GHz and 22 GHz as 
a« = log(5Q/Sjc)/log(43/22), where S Q and S K 
are flux of components at 43 GHz and 22 GHz. 
As shown in Figure El the spectrum of CI is flat 
( a Q rv> 0) while those of C2 and C3 are steep 
(o;Q ~ —0.9). It seems natural to interpret CI 
as the radio core and C2 and C3 as optically-thin 
jet components. Hereafter, we regard CI as the 
reference position. 

Figure |4] shows the change of the peak posi- 
tion of C3. The motion of C3 is mostly in the 
north-south direction, but also shows the motion 
in the east-west direction. To describe the posi- 
tional change of C3, firstly we set a coordinate as 
shown in Figure [5] We define x\\ axis to be parallel 
to the line between the position of CI and C3 in 
2003 November 20 (position angle = 161.4°) and 
x± axis to be perpendicular to xu. Then, we fit- 
ted the positional change of C3 using a quadratic 
function in consideration for the change in speed 
for each axis (Figure [5|). Position of C3 in xn and 
x±_ coordinates is parameterized as, 

X|l (t) =at 2 +bt + c, (4) 
x±(t) =dt 2 + et + /, (5) 

where the units of X|| and x± are in mas, t is the 
time from 2003 November 1 in years. The best fit 
parameters are, a = (0.028 ± 0.006), b = (0.076 ± 
0.030), c = (0.147 ± 0.046), d = (0.023 ± 0.005), 
e = (0.044 ± 0.028), and / = (0.005 ± 0.039). The 
apparent speed and acceleration can be calculated 
from these parameters. The apparent speed of C3 
is changing from 0.09±0.04 mas yr -1 to 0.41±0.07 
mas yr" 1 (or /3 app = 0.10±0.05 to 0.47±0.08) be- 
tween 2003 November 20 and 2008 November 27 
with acceleration rate of 0.07 ± 0.01 mas yr -2 (or 
^a PP = 0.08 ± 0.02 yr- 1 ). The averaged /3 app of 
C3 from 2003 November 20 to 2007 November 2, 
when C3 is not identified in Paper I, is derived to 
be (0.23±0.06)c. Although Paper I discussed pos- 
sible detection of relativistic flow with /3 app ~ 2.3 
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during the middle of 2007, we cannot find such a 
fast motion at earlier epochs from high-resolution 
VLBA data. Thus, we can rule out the existence 
of relativistic flow with /3 app ~ 2.3 during 2006 - 
2007. In Figure [J51 the curved line shows the con- 
straint on the jet speed and jet angle to the line 
of sight, based on the observed projected speed of 
C3 (0.10c- 0.47c), to gether with t he jet speed and 
jet angle adopted in lAbdo et all (|2009h . There is 
a clear discrepancy between the measured motion 
and predicted jet speed from both one-zone SSC 
model and the deceleration model. 

The averaged apparent speed of C3 from 2007 
November 2 to 2008 November 27 is derived as 
(0.43 ± 0.09)c from interpolation. This value is 
somewhat larger than the result of Paper I (/3 app = 

0. 23 ± 0.01 on average between 2007 October 24 
and 2009 April 24). This inconsistency might be 
caused by (1) the difference of coordinate setting 
and time-dependence of fitting function, and/or 
(2) under-estimation of positional error in Paper 

1. Focusing on the epoch overlapping with Paper 
I (2007-2008), the separation of C3 is systemat- 
ically larger than that in Paper I by ~ 0.1 mas 
as if expected from the position shift of reference 
position (peak position of CI) due to the optical 
depth. 

C2 is too faint to accurately determine its 
peak position, and higher dynamic range data are 
needed to identify and discuss it's kinematics. 

4.3. Lightcurves 

Figurc[7]shows the lightcurve of total CLEANed 
flux and the flux density of each component 
at 43 GHz, together with Metsahovi 37-GHz 
lightcurve. The flux of each component is defined 
as the integrated flux of Gaussian component, 
which is identified with the Gaussian models ob- 
tained by the Difmap modclfit. An example of 
the result of Gaussian fitting is shown in Fig- 
ure [3 Overall trends of the 43-GHz VLBA total 
CLEANed flux and Metsahovi 37-GHz lightcurve 
are very similar, and the flux density of Metsahovi 
at 37 GHz is always brighter by ~ 5 Jy than that 
of VLBA at 43 GHz throughout all epochs. This 
suggests that the flux variability essentially arises 
from the central pc region. 

From 2002 to 2008, the lightcurve of total 
VLBA flux density shows moderate increase. Af- 



ter 2008, rapid flux increase is clearly seen. C3 
shows a significant flux density increase, which is 
consistent with Paper I. CI also shows a flux den- 
sity increase, although not obvious in Paper I. This 
could be because of large optical depth at 22 GHz. 
The increasing tendencies of CI and C3 flux den- 
sity resemble that of total CLEANed flux. This 
indicates that the radio flare in central ~ 1 pc 
region originates in CI and C3, invoking the pos- 
sibility that the GeV 7-ray emission originates in 
these innermost components. This view is con- 
sistent with the 7-ray em itting zone derived by 
Leon-Tavares et al. ( 201lh (distance from the core 
< lpc). 

5. SED MODELING AND INTERPRE- 
TATION 

In this section, we try to reproduce the broad- 
band flux from 7-ray emitting region of 3C 84 
based on the quantities of C3 obtained as in 2J In 
Figure [51 we show the observed SED of 3C 84 and 
one-zone SSC model fit to the observed SED using 
the size, flux, and /3 app of C3 component measured 
by the VLBA obse rvations. Th e flux data are the 
same as the ones in lAbdo et all (|2009l ) but we add 
the 43-GHz flux of C3 measured on 2008 August 
27. We assume that 7-ray and VLBA 43-GHz flux 
are emitted from the same sub-pc region. Optical 
data (afte r the extract i on of host galaxy contribu- 
tion; c.f., lAbdo et al. ( 20091 )) is also used as flux 
upper limit of C3. The /3 app at 2008 August is 
0.44 (in Figure [6j . Assuming the jet viewing an- 
gle as 25° (jAbdo et all (j2009h . the intrinsic jet 
speed (3 can be estimated as 0.54, which corre- 
sponds to r = 1.19 (S = 1.65). The Gaussian- 
fitted model size of C3 is ~ 0.2-0.3 mas. To fit 
the 43-GHz flux, the size of 7-ray emitting region 
is adopted as 1.45 x 10 17 cm (0.133 mas) which 
is slightly smaller than the C3 model size. The 
detailed numerical treatment of synchrotron and 
inverse Compton scattering processes is shown in 
Kino et al.l ([2002J) and references therein. 

If we adopt 7 m in ~ 10 3 as is the case in 



_2inin 

m, ■ 



Abdo et al. (20091), 43 GHz lies below the fre- 



quency of low energy cutoff and therefore the spec- 
tral index should be a® =1/3 (vf v oc 4/3). This 
is inconsistent with the observed spectral index be- 
tween 22 and 43 GHz a.% = -0.87 ± 0.3 as shown 
in Figure [3] Then, we adopt 7 m ; n = 1 to raise the 
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cutoff frequency. This fit is indicated by a solid 
line in Figure [9l The model spectrum requires an 
extremely hard electron energy inde x s = 1.2 occ a- 



sional ly seen in TeV blazars (e.g., iKataoka et al 
( 2000f )) if we try to reproduce the flux data at 
other wavelengths. The synchrotron opacity be- 
comes unity at ~ 30 GHz, and the resultant model 
spectrum shows q« ~ —0.05, which is still incon- 
sistent with the observed spectral index. Hence 
it seems difficult to attribute the broadband SED 
to the one-zone SSC emission from C3. This indi- 
cates that some other components CI and/or C2 
also should be taken in account for the candidates 
for 7-ray emitting region. 

6. SUMMARY 

Following the detection of 3C 84 by Fermi /L AT 
in 2008, we discussed the kinematics of the newly 
formed jet component C3 using our VERA radio 
observations for the period 2006 June - 2009 April 
in Paper I. In this paper, we further explored the 
kinematics and lightcurves of CI, C2, and C3 us- 
ing archival 43-GHz VLBA data spanning a time 
interval from 2002 January to 2008 November, i.e. 
before the VERA observations. Summary and dis- 
cussions are as follows: 

1. In the multi-epoch high resolution images 
obtained by VLBA at 43 GHz (Figure [2), 
we find that C3 has been visible since 2003 
November. We succeed in obtaining a strong 
constraint on the time of the C3 birth before 
2003 November. 

2. The lightcurvc of VLBA at 43 GHz shows 
a trend similar to that of Metsahovi at 
37 GHz, indicating that the increase of ra- 
dio flux density arises from the region within 
the central ~ 1 pc. In particular, the rapid 
increase of flux density has been seen since 
the middle of 2008. Both CI and C3 show 
the largest change of flux density for this 
flare. This indicates that the radio flare in 
the central ~ 1 pc region originates in both 
CI and C3. The flux density increase of 
CI was not prominent in Paper I, possibly 
because of large optical depth at 22 GHz. 

3. The apparent speed of C3 with respect to 
CI changes from (0.10 ± 0.05)c to (0.47 ± 
0.08)c between 2003 November 20 and 2008 



November 27 with acceleration rate /3 app = 
0.08 ± 0.02 yr _1 . No superluminal motion is 
detected with VLBA observations at 43 GHz 
before 2008 November. 

4. The one-zone SSC model using measured ap- 
parent speed of C3 is fitted to the observed 
broadband spectrum, in particular, the 43- 
GHz and GeV 7-ray fluxes. Our model fit is 
difficult to reproduce the optically-thin radio 
spectrum ~ —0.9 that is measured by 
the VLBA observations. This indicates that 
it seems difficult to attribute the broadband 
SED to the one-zone SSC emission from C3. 
Therefore, the other components should not 
be excluded as candidates for 7-ray emitting 
regions. 
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Fig. 1. — Estimation of position measurement 
accuracy for C3. Each point represents the differ- 
ence of C3 position measured between two epochs 
closely separated in time (within 10 days), which is 
normalized by the beam size averaged between the 
pair. Thick bars correspond to the standard de- 
viation of the scatter for right ascension and dec- 
lination direction, and arc estimated to be 0.255 
6>beam and 0.166 6> bc am, respectively. 
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Fig. 2.— The VLBA 43-GHz images from 2002 
January 7 to 2008 November 27. Image sizes are 4 
x 4 mas 2 (1 mas corresponds to 0.353 pc). All im- 
ages are convolved with the same restoring beam 
of 0.3 mas in diameter. Contour levels are plot- 
ted at the level of 68.4 x n 2 (n = 1,2,3,4,) mJy 
beam -1 . The lowest contour is the highest image 
noise r.m.s. among all images. 



10.0 




20 50 100 

Frequency [GHz] 

Fig. 3. — Spectral indices of CI, C2 and C3 be- 
tween 22 GHz and 43 GHz at 2008 November 27. 
We assume that flux error is 10% of the flux den- 
sity of each component. Spectral indices of CI, 
C2, and C3 are 0.06 ± 0.30, -0.92 ± 0.30, and 
-0.87 ± 0.30, respectively. 
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Fig. 4. — Peak position of C3 for all epochs, su- 
perposed on the contours of 43-GHz intensity dis- 
tribution on 2008 November 27. The dashed line 
traces the average positional change of C3. We 
set a; || and x± coordinates to evaluate the motion 
parallel and perpendicular to the initial position 
angle of C3, 161.4° from CI. 
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Fig. 5. — The change of relative angular distance 
between CI and C3 in x\\ and x± coordinates as a 
function of time from 2003 November 1 in years. 
Red and blue points show and x± for each 
epoch. Orange and cyan lines are fitted functions. 
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Fig. 6. — Constraint on the bulk Lorentz factor T 
and the jet angle to the line of sight 9. The gray- 
filled area (corresponding 0.10 < /? app < 0.47) is 
allowed for the apparent speed range of C3 be- 
tween 2003 November and 2008 November. The 
white circle and gray-filled circle show the best-fit 
parameters for one-zone SSC model (r — 1.8 with 
9 = 25°) and for deceleration jet model (as shown 
by the arrow, jet decelerates from T max — 10.0 to 
r m in = 2.0 over a dist ance of 5 x 10 17 cm (0.16 
pc), with 9 = 20°) from lAbdo etaD (|2009h . 
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Fig. 7. — VLBI and single dish lightcurves 
of 3C 84. Pink and magenta points indi- 
cate the flux density of Metsahovi at 37 GHz 
(jLeon-Tavares et al.|[20Tlh and VLBA at 43 GHz, 
respectively. Flux densities of CI, C2, and C3 are 
indicated by blue, green, and red points, respec- 
tively. 
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Fig. 8. — An example of model fitting result for 
VLBA 43-GHz image. Black contour shows the 
brightness distribution. Blue, green, and red sym- 
bols correspond to CI, C2, and C3, respectively. 
Black symbols are other minor components. 
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Fig. 9.— Overall SED of NGC 1275 con- 
structed with mult iband data (black and magenta; 
Abdo et al.l (|2009t )) and the VLBA 43-GHz flux of 
C3 (red square; this work). The SED is fitted 
to VLBA 43-GHz flux of C3 (red square), opti- 
cal/UV (Swift/UVOT; magenta open circle), and 
7-ray (Fermi; magenta point) data with a one-zone 
SSC model. We obtained the model parameters 
as R = 1.45 x 10 17 cm (0.133mas), B = 0.37 G, 



1.2, 



In 



1. 



q e = 5 x 10 

7max = 7.5 x 10 4 , and S = 1.65, where R is the 
comoving radius of the jet emitting region, B is 
the mean magnetic field in the radiating plasma, 
<Ze, s, 7 m i n , 7 max are the normalization factor, the 
power-law index, maximum an d minimum Lorent z 
factor for injection spectrum in lKino et al.l ([20021 ). 
Doppler factor S is derived from the measured 
value of /3 app = 0.44 for C3 at 2008 August and 
the assumed viewing angle 8 = 25 deg (yielding T 
= 1.19). The flux measured at 43 GHz on 2008 
August 27th is 2.67 ± 0.32 Jy. 
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Table 1 
Observational Information. 



Arr£iy 


"Fin orb 


Band a 


Tirnp On 

Sonrcp (min) 


VLBA 


2002 Jan 7 


Q 


32 


VLBA 


2002 Apr 5 


Q 


704 


VLBA 


2003 Jan 12 


Q 


25 


VLBA 


2003 Nov 20 


Q 


6 


VLBA 


2004 Oct 2 


Q 


6 


VLBA 


2005 Jan 25 


Q 


7 


VLBA 


2005 Oct 4 


Q 


6 


VLBA 


2006 Jan 7 


Q 


9 


VLBA 


2006 Mar 22 


Q 


9 


VLBA 


2006 Apr 7 


Q 


6 


VLBA 


2006 May 08 


Q 


6 


VLBA 


2006 May 26 


Q 


6 


VLBA 


2006 Jul 26 


Q 


9 


VLBA 


2006 Sep 1 


Q 


6 


VLBA 


2006 Sep 22 


Q 


6 


VLBA 


2007 Jan 4 


Q 


9 


VLBA 


2007 May 24 


Q 


9 


VLBA 


2007 Jun 15 


Q 


6 


VLBA 


2007 Aug 7 


Q 


5 


VLBA 


2007 Nov 2 


Q 


G 


VLBA 


2008 Apr 3 


Q 


6 


VLBA 


2008 Apr 16 


Q 


15 


VLBA 


2008 May 8 


Q 


5 


VLBA 


2008 May 23 


Q 


6 


VLBA 


2008 Jun 31 


Q 


6 


VLBA 


2008 Aug 27 


Q 


6 


VLBA 


2008 Nov 1 


Q 


6 


VLBA 


2008 Nov 27 


K,Q 


6 



a Band K = 22 GHz, Q = 43 GHz. 
b Total on source time in min. 
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Table 2 
Dispersions of C3 positions. 



Pairs a 


nmean b 


A RA c 


A Dec c 




(mas) 


(beam) 


(beam) 


2004 Oct 2 / 2004 Oct 3 


0.31 


-0.01 


-0.01 


2005 Oct 4 / 2005 Oct 12 


0.33 


-0.20 


0.09 


2005 Oct 12 / 2005 Oct 18 


0.33 


0.16 


0.06 


2006 Jan 6 / 2006 Jan 7 


0.27 


-0.04 


-0.03 


2006 Mar 22 / 2006 Mar 23 


0.27 


0.07 


-0.11 


2006 Apr 7 / 2006 Apr 15 


0.32 


0.26 


0.07 


2006 May 8 / 2006 May 12 


0.34 


0.06 


-0.19 


2006 Sep 22 / 2006 Sep 28 


0.28 


-0.01 


0.05 


2006 Dec 28 / 2007 Jan 4 


0.30 


0.09 


0.01 


2007 Aug 7 / 2007 Aug 17 


0.33 


-0.22 


-0.14 


2008 Nov 1 / 2008 Nov 11 


0.36 


-0.17 


-0.14 



a Pairs of near epochs with their separation < 10 days (for 
which the motions of components are negligible) 

b Bcam sizes averaged for the major - minor axis and two 
epochs. 

c Diffcrences of relative right ascension and declination po- 
sition of C3 between two epochs. 
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Table 3 

Epoch, rms, synthesized beam size, peak brightness, total CLEANed flux for all image. 



Epoch 


rms 


Beam a 


T b 

-tpeak 


S c 




(mjy beam -1 ) 


(mas x mas, degree) 


(Jy beam -1 ) 


(Jy) 


9009 Tan 7 


8 ^ 
o.o 


n qi y n 1Q on 


1 07 


Q AO _|_ fl QC 

0.4:0 zn u.oo 


2002 Anr 5 


6 6 


25 x 14 -33 7 


U. £7U 


• >._■/ _l_ u.uu 


2003 Tan 12 


10.6 


39 x 16 3 1 


1.00 


4.07 ± 0.41 


2003 Nov 20 


23 7 


0^3 x 18 33 6 


1 1 7 

1 . 1 1 




2004 Ort 2 


1 3 


49 x 14 -23 1 


1 59 

1- .OtJ 


A 49 + A A 

4: . 4: Z; _l_ U . 11 


2005 Tan 25 

ZjUWU J dill ZiU 


12.8 


29 x 17 -22 2 


1.04 


a 1 n + n 41 


2005 Ort 4 


1 8 8 


1 7 x 46 -34 4 


1 49 


4 (if) + n 4« 

4:.UU —1— U.t:U 


200fi Tan 7 

ZjVJUU tJ till 1 


39 8 


27 x 1 9 -5 7 


1.44 


4 4fl _|_ n 44 

t: . t: U 1 U . 11 


900fi A/Tnr 99 


1 7 3 


Of) y ^ fi 8 

U..ZU A U.OO > U.O 


1 1 8 

J — LO 


q 79 _l n 07 

O. ( Z ZC U.O ( 


900fi Anr 7 


uo.i 


n w y niR 14 4 


1 fiO 


A 97 _L A AO 


900fi 1Y4W & 

ZiUUU 1V1 ay O 


4^ 3 


U.ZU A U.^O , Z,O.U 


1 7^ 

1 . ( o 


c f)7 _l n 

ej . U I ZC U.U 1 


2006 May 26 


13.4 


0.48 x 0.17 , -33.0 


1.30 


4.34 ± 0.43 


2006 Jul 26 


36.3 


0.44 x 0.23 , -4.2 


1.58 


4.20 ± 0.42 


2006 Sep 1 


15.5 


0.16 x 0.48 , -34.2 


1.32 


3.77 ± 0.38 


2006 Sep 22 


34.2 


0.25 x 0.19 , -13.0 


1.45 


4.94 ± 0.49 


2007 Jan 4 


38.7 


0.28 x 0.19 , -15.8 


1.34 


4.48 ± 0.45 


2007 May 24 


20.6 


0.22 x 0.19 , -11.7 


0.99 


3.93 ± 0.39 


2007 Jun 15 


39.4 


0.26 x 0.55 , 14.5 


2.04 


4.37 ± 0.44 


2007 Aug 7 


45.0 


0.51 x 0.23 , 16.9 


1.73 


5.43 ± 0.54 


2007 Nov 2 


39.2 


0.18 x 0.53 , -1.9 


2.16 


4.62 ± 0.46 


2008 Apr 3 


33.7 


0.14 x 0.47 , -22.5 


1.85 


5.82 ± 0.58 


2008 Apr 16 


14.4 


0.45 x 0.19 , -15.3 


1.93 


4.93 ± 0.49 


2008 May 8 


41.0 


0.49 x 0.14 , -11.1 


2.00 


5.97 ± 0.60 


2008 May 23 


37.8 


0.46 x 0.15 , 20.9 


1.43 


4.36 ± 0.44 


2008 Jul 31 


46.5 


0.51 x 0.15 , -12.2 


2.02 


6.33 ± 0.63 


2008 Aug 27 


34.3 


0.54 x 0.20 , 27.3 


2.46 


6.89 ± 0.69 


2008 Nov 1 


57.1 


0.53 x 0.18 , -1.0 


2.85 


9.19 ± 0.92 


2008 Nov 27 


61.3 


0.45 x 0.14 , 7.6 


2.22 


8.90 ± 0.89 



a Major axis, minor axis, and position angle of synthesized beam. 
b Peak brightness for each image. 

c Total CLEANed Flux and its error for each image. We assumed that the flux calibration error is 10% of 
flux density. 
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Table 4 

model type and flux of fitted ci component. 



Epoch 


ModeP 


C b 
^mod 






(Jy) 



2002 Jan 7 2002 


c 


2 


.13 


± 


0. 


22 


2002 Apr 5 2002 


e 


1 


.16 


± 


0. 


.12 


2003 Jan 12 2003 


e 


2 


.04 


± 


0. 


.21 


2003 Nov 20 2003 


e 


2 


.68 


± 


0. 


.29 


2004 Oct 2 2004 


e 


2 


.15 


± 


0. 


.22 


2005 Jan 25 2005 


c 


2 


.27 


± 


0. 


24 


2005 Oct 4 2005 


c 


1 


.93 


± 





.22 


2006 Jan 7 2006 


c 


2 


.13 


± 


0. 


22 


2006 Mar 22 2006 


e 


1 


.59 


± 





,17 


2006 Apr 7 2006 


e 


2 


.61 


± 





.30 


2006 May 8 2006 


c 


2 


.09 


± 





22 


2006 May 26 2006 


c 


1 


.34 


± 


0. 


.13 


2006 Jul 26 2006 


c 


1 


.64 


± 


0. 


17 


2006 Sep 1 2006 


e 


1 


.74 


± 





18 


2006 Sep 22 2006 


c 


1 


.77 


± 





,18 


2007 Jan 4 2007 


c 


2 


.01 


± 





22 


2007 May 24 2007 


e 


1 


.69 


± 


0. 


,18 


2007 Jun 15 2007 


c 


1 


.97 


± 


0. 


,20 


2007 Aug 7 2007 


c 


1 


.72 


± 


0. 


,18 


2007 Nov 2 2007 


c 


2 


.62 


± 


0. 


.28 


2008 Apr 3 2008 


c 


2 


.22 


± 





23 


2008 Apr 16 2008 


c 


1 


.89 


± 





,19 


2008 May 8 2008 


c 


2 


.46 


± 


0. 


.25 


2008 May 23 2008 


e 


2 


.15 


± 


0. 


.23 


2008 Jan 31 2008 


c 


2 


.67 


± 


0. 


.28 


2008 Aug 27 2008 


e 


3 


.37 


± 


0. 


35 


2008 Nov 1 2008 


e 


4 


.36 


± 





,47 


2008 Nov 27 2008 


c 


3 


.49 


± 





37 



a Shape of component models, 'e' is ellipti- 
cal Gaussian, V is circle Gaussian function. 

b Flux and its error for models. Error is 
estimated as the addition of calibration error 
(typically 10% of component) and image rms 
of each epoch. 



14 



Table 5 

Relative position from CI, model type, and flux of fitted C2 component. 



Epoch 


Relative R.A. a 
(mas) 


Relative Dec. a 
(mas) 


Model b 


(Jy) 


9009 Tan 7 9009 


^ 
-u.oo 


78 
-u. t o 


c 


1 . _L y ZC U . 14: 


9009 Am K 9009 
ZUUZ r\.pi ZUUZ 


ftft 

-u.uu 


88 
-U.oo 


c 




9001 Tan 1 9 9001 
ZUUO J cLIl ±Z ZUUO 


ft^ 

-u.uo 


77 
-U. i t 


e 


u.oo zlz u.iu 


9001 TVrnr 90 9001 
ZUUO 1NOV ZU ZUUO 


1ft 

-u.ou 


7ft 
-U. ( u 


c 


n 4- n no 
u.oy ziz u.uy 


900/1 Ont 9 900/1 
ZUU^ UlL z ZUU^ 


70 
-U. ( u 


01 

-u.yo 


e 


u.yo zlz u.io 


900^ Tan 9^ 900^ 
ZUUO J ail ZO ZUUO 


81 
-U.Ol 


01 

-u.yo 


c 


U.^Z ZIZ u.uo 


900^ DpT A 900^ 

ZUUO WLL ^± ZUUO 


^1 
-u.oo 


Q^i 
-u.yo 


c 


i qo _i_ n i c 

1 .0/i ZLZ U . 1 <J 


900fi Tan 7 900fi 
ZUUU Jail / zuuu 


R9 
-u.uz 


1 1 9 


e 


U.OH: zlz u.uo 


900fi Mar 99 900fi 

ZUUU IvlcLL ZZ ZUUU 


ft/L 


111 

-1 . lO 


c 


U.Ol ZIZ U. 1U 


900R Am 7 900ft 
ZUUU r\pi 1 ZUUU 


89 
-u.oz 


1 9^> 

- _L . ZO 


e 


70 -I- 1 1 

U.IU ZLZ U. 1 1 


900ft A/fair 8 900ft 

zuuu lviay o zuuu 


^7 
-U.Ol 


11ft 
-1 . 1U 


c 


1 94 -I- 1 f\ 

1 .Z^ ZLZ U. lO 


2006 May 26 2006 


-0.53 


-1.11 


e 


1.44 ± 0.18 


2006 Jul 26 2006 


-0.73 


-1.07 


e 


0.94 ± 0.14 


2006 Sep 1 2006 


-0.55 


-1.15 


e 


1.19 ± 0.16 


2006 Sep 22 2006 


-0.56 


-1.13 


e 


1.35 ± 0.20 


2007 Jan 4 2007 


-0.65 


-1.19 


e 


0.66 ± 0.11 


2007 May 24 2007 


-0.57 


-1.40 


e 


0.43 ± 0.06 


2007 Jun 15 2007 


-0.74 


-1.30 


e 


0.61 ± 0.08 


2007 Aug 7 2007 


-0.96 


-1.27 


c 


0.63 ± 0.06 


2007 Nov 2 2007 


-0.69 


-1.34 


c 


0.61 ± 0.10 


2008 Apr 3 2008 


-0.63 


-1.24 


e 


0.82 ± 0.20 


2008 Apr 16 2008 


-0.57 


-1.22 


c 


0.82 ± 0.13 


2008 May 8 2008 


-0.74 


-1.11 


c 


0.53 ± 0.16 


2008 May 23 2008 


-0.60 


-1.16 


c 


0.57 ± 0.13 


2008 Jan 31 2008 


-0.83 


-1.07 


c 


0.35 ± 0.04 


2008 Aug 27 2008 


-0.35 


-1.24 


d 


0.51 ± 0.05 


2008 Nov 1 2008 


-0.74 


-1.49 


d 


0.37 ± 0.04 


2008 Nov 27 2008 


-0.71 


-1.47 


d 


0.36 ± 0.04 



a Rclativc R.A. and Dec. angular distances between CI and C2. 

b Shape of component models, 'c' is elliptical Gaussian, 'c' is circle Gaussian, 
and 'd' is delta function. 

c Flux and its error for models. Error is estimated as the addition of calibra- 
tion error (typically 10% of component) and image rms of each epoch. 
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Table 6 

Relative position from CI, model type, size, and flux of fitted C3 component. 



Epoch 


Relative R.A. a 


Relative Dec. a 


Model b maj x min , P.A. c 


c d 




(mas) 


(mas) 


(mas x mas, deg) 


(Jy) 



2002 Jan 7 2002 

2002 Apr 5 2002 

2003 Jan 12 2003 



iUUO 1\UV ZU iUUO 


n 
u 


(U 


n 
-u 


1 9 
1 z 


c 


U.IO A u.io , u.u 


n 
u 


uo 


4- 


(1 
u 


(17 
u t 




n 
u 


no 

U -J 


-u 


zo 


c 


U.Zi A U.UU , uu.o 


1 

1 


1 A 


4- 


(1 


1 1 

1 1 


9(10^. Tan 9^ 900^. 
zuuo j an zo zuuo 


n 
u 


1 9 
1Z 


n 
-u 


9 J. 

Z4: 


c 


U.UO A u.uo , u.u 


n 
u 


oz 


4- 


n 
u 


uo 


900S Oct 4 900S 


n 

u 


1 '-i 


-0 


4fi 




93 x 09 97 9 


I 


34 


-j- 




u 


14 


2006 Jan 7 2006 





14 


-0 


47 


c 


0.13 x 0.13 , 0.0 





67 


± 





07 


2006 Mar 22 2006 





13 


-0 


50 


c 


0.23 x 0.23 , 0.0 


1 


01 


± 





11 


2006 Apr 7 2006 





14 


-0 


48 


c 


0.06 x 0.06 , 0.0 





81 


± 





08 


2006 May 8 2006 





13 


-0 


50 


c 


0.23 x 0.23 , 0.0 


1 


72 


± 





19 


2006 May 26 2006 





12 


-0 


52 


e 


0.43 x 0.23 , 1.9 


1 


71 


± 





19 


2006 Jul 26 2006 





03 


-0 


45 


c 


0.20 x 0.20 , 0.0 


1 


01 


± 





12 


2006 Sep 1 2006 





12 


-0 


54 


c 


0.32 x 0.32 , 0.0 


1 


18 


± 





14 


2006 Sep 22 2006 





11 


-0 


55 


c 


0.29 x 0.29 , 0.0 


1 


58 


± 





19 


2007 Jan 4 2007 





04 


-0 


59 


e 


0.45 x 0.21 , -61.5 


1 


37 


± 





21 


2007 May 24 2007 





15 


-0 


84 


c 


0.22 x 0.22 , 0.0 





88 


± 





11 


2007 Jun 15 2007 





17 


-0 


80 


e 


0.00 x 0.00 , 0.0 





91 


± 





09 


2007 Aug 7 2007 





13 


-0 


95 


e 


0.55 x 0.29 , 17.7 


1 


89 


± 





25 


2007 Nov 2 2007 





06 


-0 


98 


c 


0.43 x 0.28 , 45.0 


1 


74 


± 





27 


2008 Apr 3 2008 





05 


-1 


00 


e 


0.34 x 0.27 , -43.1 


1 


55 


± 





20 


2008 Apr 16 2008 





10 


-1 


09 


e 


0.44 x 0.26 , 25.7 


1 


62 


± 





18 


2008 May 8 2008 





09 


-1 


08 


e 


0.34 x 0.15 , 68.1 


1 


61 


± 





19 


2008 May 23 2008 





07 


-1 


07 


e 


0.58 x 0.33 , 16.2 


2 


04 


± 





31 


2008 Jan 31 2008 





06 


-1 


11 


e 


0.72 x 0.34 , 47.0 


2 


91 


it 





44 


2008 Aug 27 2008 





07 


-1 


11 


e 


0.65 x 0.27 , 8.5 


2 


67 


± 





32 


2008 Nov 1 2008 





05 


-1 


23 


e 


0.58 x 0.44 , 45.0 


3 


76 


± 





53 


2008 Nov 27 2008 





08 


-1 


27 


e 


0.67 x 0.34 , 43.1 


3 


57 


± 





58 



a Relative R.A. and Dec. angular distances between CI and C3. 

b Shape of component models, 'e' is elliptical Gaussian, V is circle Gaussian function. 

c Major axis, minor axis, and position angle of models. (? ma j and 9 m - m are not defined for delta 
models. P.A. is not defined for circle and delta models. 

d Flux and its error for models. Error is estimated as the addition of calibration error (typically 
10% of component) and image rms of each epoch. 
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